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Abstract — Distributed opportunistic sclieduling is studied for 
wireless ad-hoc networks, where many links contend for one 
channel using random access. In such networks, distributed op- 
portunistic scheduling (DOS) involves a process of joint channel 
probing and distributed scheduling. It has been shown that under 
perfect channel estimation, the optimal DOS for maximizing the 
network throughput is a pure threshold policy. In this paper, this 
formalism is generalized to explore DOS under noisy channel 
estimation, where the transmission rate needs to be backed off 
from the estimated rate to reduce the outage. It is shown that the 
optimal scheduling policy remains to be threshold-based, and that 
the rate threshold turns out to be a function of the variance of the 
estimation error and be a functional of the backoff rate function. 
Since the optimal backoff rate is intractable, a suboptimal linear 
backoff scheme that backs off the estimated signal-to-noise ratio 
(SNR) and hence the rate is proposed. The corresponding optimal 
backoff ratio and rate threshold can be obtained via an iterative 
algorithm. Finally, simulation results are provided to illustrate the 
tradeoff caused by increasing training time to improve channel 
estimation at the cost of probing efficiency. 

I. Introduction 

Channel-aware scheduling for achieving rich diversities 
inherent in wireless communications has recently emerged as 
a promising technique for improving spectral efficiency in 
wireless networks. Most existing studies along this line require 
centralized scheduling (see [5] and the references therein), and 
httle work has been done on developing distributed algorithms 
to harvest diversity gains for ad hoc communications. Unlike 
centralized cases, in ad hoc communication networks, each 
link has no knowledge of other links' channel conditions, mak- 
ing it very challenging to carry out channel-aware distributed 
scheduling. 

Some initial steps have been taken by several of the authors 
to develop distributed opportunistic scheduling (DOS) to reap 
multiuser diversity and time diversity in wireless ad hoc 
networks [7]. In DOS, a successful link proceeds to data 
transmission only if the observed channel condition is "good"; 
otherwise, it may skip the transmission, and let all the links re- 
contend for the channel. Intuitively speaking, because different 
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links in different time slots experience different channel con- 
ditions, it is likely that after further probing, the channel can 
be taken by a link with a better channel condition, resulting in 
possible higher throughput. In this way, the multiuser diversity 
across links and the time diversity across slots can be exploited 
in a joint manner 

Despite the insightful analytical results reported in [7], its 
key assumption is that perfect channel state information (CSI) 
is known at the receiver/transmitter However, in practical 
scenarios, channel conditions are often estimated using noisy 
observations. Therefore, it is of great interest to study DOS 
under noisy channel estimation. In centralized scheduling 
schemes, it has been shown that the estimated signal-to-noise 
ratio (SNR) is always larger than the "actual SNR" [6]. Thus, 
if the data were transmitted using the estimated rate, there 
would always be an outage. To reduce the outage probability, 
the transmission rate has to back off from the estimated rate. 
Therefore, unlike the perfect estimation case, the optimal 
scheduling policy hinges on the backoff rate. 

In this work, we generalize [7] to the scenario with imper- 
fect channel estimation and show that the optimal schedul- 
ing policy for the noisy channel estimation case remains a 
threshold structure. However, the threshold turns out to be a 
function of the variance of the channel estimation error, and 
furthermore, it is a functional of the backoff rate function. 
Since the optimal backoff rate function is difficult to obtain, 
we propose a suboptimal linear backoff scheme. We show that 
the corresponding optimal backoff ratio and rate threshold can 
be obtained via an iterative numerical algorithm. Simulation 
results are provided to show that DOS achieves significant 
throughput gain in the presence of noisy channel estimation, 
especially in the low SNR region. 

II. System Model and Background 

We consider a single-hop ad hoc network with M links, 
where link m contends for the channel with probability p„i, 
m — 1 , . . . , M. A collision model is assumed for random 
access, where channel contention of a link is said to be suc- 
cessful if no other links transmit at the same time. Accordingly, 
the overall successful contention probability, ps, is then given 
by Z]m=i(Pm rii^^mll ~ P^)) ■ ^Icar that the number of 
slots (denoted as K) for a successful channel contention is a 
Geometric random variable (rv.), i.e., K ~ Geometric{ps). 
Let T denote the duration of mini-slot for channel contention, 
and T the data transmission time. It follows that the random 
duration corresponding to one round of successful channel 
contention is Kt, with expectation r/ps. 



Let s(n) denote the successful link at the n-th successful 
channel contention. The corresponding received signal is given 
by: 

Ys{n) (n) = ^/phs(„) (n)Xs(„) (n) + {n), (1) 

where p is the normalized receiver SNR, /is(„)(n) is the 
channel coefficient for link s{n), Xg(^n)(n) is the transmitted 
signal with _E |||Xs(„)(r7,)|p} = 1 and ^is(n){n) is additive 
white noise with i.i.d. CA/^(0, 1). 

To simplify the exposition, we consider a homogeneous 
network in which all links have the same channel statistics, 
and are subject to Rayleigh fading, i.e., hs(n){n) follows a 
complex Gaussian distribution CA/'(0, 1). In what follows, we 
drop the subscripts to simplify the notation and use hn to stand 
for hg(n-f{n) where it is clear from the context. Similarly, we 
use Yn, Xn and ^„ to denote (n), Xs(„) (n) and /x^(„) (n). 

We consider the continuous rate case, assuming that the 
instantaneous rate is given by the Shannon channel capacity, 
i.e., 

Rn = log(l + p\hn\^) nats/s/Hz, 

provided that the channel can be perfectly estimated. 

In [7], we have studied DOS with perfect CSI. Specifically, 
we have shown that the problem can be cast as a maximal 
rate of return problem in optimal stopping theory [3], where 
the rate of return is the average network throughput, x, and is 
determined by the stopping time N: 

E[RnT] 



E[Tn] 



(2) 



where T„ = J2]=i + ^ '^^e total system time including 
the contention time and the data transmission time. Note that 
iV is a stopping time if {N — n} is JF„ -measurable, where 
!Fn is the cr-field generated hy {{p\hj\'^ , Kj), j = 1,2, ... ,n}. 

We show that the optimal DOS maximizing the throughput 
is given by the optimal stopping rule, N*, that solves the 
maximal rate of return problem in (2), i.e.. 



N 



» a EIRnT] 

— arg max ' . 
NeQ E[Tn\ 



E[RnT] 



where 



Q^{N ■N>1, E[Tn] < oo}. 



(3) 



(4) 



It has been shown in [7] that the optimal stopping rule N* 
for DOS exists, and is given by 



N* = min{n > 1 : i?„ > x*}. 



(5) 



Furthermore, the maximal throughput x* is an optimal thresh- 
old, and is the unique solution to 

XT 



E{R-x)+ = —, 
PsT 



(6) 



where i? is a r.v. having the same distribution as Rn- 

For example, if Rn = log(l+p|/in|^) and /i„ has a complex 
Gaussian distribution CAf{0, 1), it can be shown that 



'1\ /exp(ar^\p, 
. =exp(-ji.,^-^j-, 



(7) 



where S = t/T, and Ei{x) is the exponential integral function 
defined as Ei(x) 4 /~ ^-^S^dt. 

III. DOS Under Noisy Channel Estimation 



Needless to say, in practical systems, hn has to be estimated 
using training signals (e.g. embedded in the RTS packets). Let 
hn denote the estimation of the channel coefficient, and hn the 
estimation error. It follows that 



hn — ~t~ 



(8) 



where hn and hn are zero-mean complex Gaussian random 
variables. Suppose that the channel is estimated using a mini- 
mum mean square error (MMSE)-based estimator. It follows, 
by the orthogonality principle, that 



E[\hn\'']=E[\hn\'']+E[\hn\'']. 



(9) 



Let (3 denote the variance of the estimation error. From (9), 
we have that 



£[|/i„|2]=/3 and ^[|/j„|2] = l-/3. 



(10) 



Treating the estimation error as noise, the actual SNR at the 
receiver can be computed by [6] 



An 



■ 



1+P\hn? 



(11) 



We note that the numerator of (11), p|/inp, is the esti- 
mated SNR. Therefore, in contrast to the perfect CSI case 
where the sequence {p\hn\'^ ,n = 1,2,...} is used for dis- 
tributed scheduUng, in the noisy estimation case, n = 
1,2,.. .} serves as the basis for distributed scheduling. 
Following [6], and can be normalized as 



An 



\hnl_ 
1-p 



and Zn 



\hn[ 



(12) 



Note that both A,i and Zn have the exponential distribution 
with unit variance. Furthermore, A„ in (11) can be rewritten 
as 

PeffK 



An — 



■ apeffZn 



(13) 



where peff = (1 — /3)p and a = denote the "effective 
channel SNR" and "normalized error variance", respectively. 
It can be shown the distribution of An given A„ is given by [6] 



1 / An 1 
a I An Peff 



An 1 
An Peff 



(14) 



where I(-) is the indicator function. 



A. Optimal Stopping Rule under Noisy Channel Estimation 

It is clear that the actual SNR A„ is no greater than the 
estimated SNR pef/Xn- As a result, if the packet is transmitted 
at the estimated rate log(l -|- pe/f Xn), there would always 
be a channel outage. Therefore, the transmission rate has to 
back off from the estimate rate. Equivalently, we can back 



off the estimated SNR pe//A to a "nominated" SNR Ac (A). 
Accordingly, the instantaneous rate, i?„ , is given by 

i?„ = log (l + Ae(A„)) I (Ae(A„) < A„) . (15) 

Along the same line as in the perfect CSl case, for each 
given back-off rate function Xc{-), maximizing the average 
throughput reduces to solving the maximal rate of return 
problem in (2). 

Observe that there are at least two major differences be- 
tween the perfect estimation case and the noisy chaimel esti- 
mation case. First, the stopping rule A'^ is now defined over the 
cr-field (instead of generated by {{p\hj\'^,Kj),j = 
1,2, ...,n}. Second, the instantaneous rate, defined in 
(15), is now a random variable, and is not perfectly known 
at time n. However, it can be shown that the structure of the 
optimal scheduUng strategy remains the same, except that the 
random "reward" i?„ is replaced with its conditional expec- 
tation, Rn = E[Rn\Pn\ [3, Page 1.3] [1]. More specifically, 
define 

Q' ^{N>l:{N = n}e T'^, ^[T^r] < oo}. (16) 

Analogously, define 

Q" ^ {AT > 1 : {AT = n} e T'l, E\Tn\ < oo}, (17) 

where T'l^ is the a-field generated by {{Rj,Kj),j = 
1,2, ... ,n}. We have the following proposition. 
Proposition 3.1: 

E[RnT] _ E[RnT] 



sup 

NeQ' 



(18) 



E[Tn] NeQ" E[Tn] 
The proof follows from Propositions 2.3, 2.4, 2.5 in [1] (with 
Xn, Zn, Wn in heu of Rn, Xn,Rn)< and the proof in [7]. Due 
to space limitation, we omit the details here. 

As a result. Proposition 3.1 indicates that the optimal 
scheduling can be based solely on given by 



log (l+Xc{Xr, 







1 — oxp < 









Xr> 



Xc{Xn) 



1 

Pejf 



where we have used the fact that P ^Ac(A„) < A„|J-"^^ = 

P (^Xc{Xn) < A„|A„j due to the independence of channel 
estimations. 

Based on the above discussion, we conclude that the optimal 
scheduling policy under noisy channel estimation is a pure 
threshold pohcy: 

Af* =min{n > 1 : ^„ > a;*}, (19) 

where the optimal threshold x* can be computed from (6), 
and hence, it is the unique solution to the following fixed 
point equation: 

x = ^{x,Xa), (20) 

where 

/^Te-Mog(l+Ac) [l-exp{- 



4>(a:,Ac)^ 



1 



_6_ 

Ps 



+ e- 



(21) 



and A' can be obtained from 



log (l + Ac(A')) 



1 — exp 



A' 



Ac(A') Peff 



X. 

(22) 



B. Optimal Backoff Rate Function 

It is clear from (20) that for a given backoff rate function 
Ac(-), there is a corresponding optimal throughput x*. There- 
fore, X* is a functional of Xc{-), denoted as a;* (Ac). We are 
interested in finding the function Ac(-) that maximizes x*(Ac), 
i.e., 

A* = argmaxa;*(Ac), (23) 
XceA 

where A is the set of the admissible functions (for example, 
A can be {Ac(A) : Ac(A) > 0,V A > 0}). 

Based on the theory of calculus of variations [4], problem 
(23) is a variational problem, and the functions A*(-) are 
called extremals. However, unhke the canonical calculus of 
variations problems, in this problem, the functional x* is not 
explicitly defined on Ac. Instead, they are connected through 
a fixed point equation. Furthermore, the integral range in (21) 
is not fixed, but is a function of Ac (cf. (22)). As a result, it 
is intractable to characterize A*. 

C. A Suboptimal Backoff Rate Function 

In what follows, we propose a suboptimal backoff rate func- 
tion, which backs off the estimated SNR by a multiplicative 
ratio a, i.e., we set 



Ac(A) = (jPeffX, 



(24) 



and < (T < 1. 

It follows from (24), (21) and (22) that the optimal through- 
put X* is the solution to 



X — ^{x, a) 
= 1 — exp 



1 



log(l + ape//A>-^' + exp (^) E, {x' + ^) 



(25) 



where 



exp 



A' = 



(26) 



It is not difficult to show that x* is a continuous and 
differentiable function of a, and hence, there exists an optimal 
backoff ratio a* such that 



a* = argmaxa;*(cr). 



(27) 



dX 



It can also be shown that a* cannot be or 1 (since the 
corresponding throughput is zero). Therefore, the optimal ratio 
a* must satisfy the first order condition '^'^ }^ ^ = 0. 



D. An Iterative Algorithm for Computing a* and x*{(j*) 

Due to the complicated structure of the fixed point equation 
(25), it is not feasible to characterize a* using the first order 
condition. In what follows, we devise an iterative algorithm 
instead using fractional optimization techniques [2]. 

Specifically, we define the following functions: 



U(a, x) 



1 — exp < — 



aPeff 



a 



|log(l + apeffX') exp(-A') + exp 



Eii\' + 



1 



^ ' exp (—A'), where A' is defined in (26). 
The iterative algorithm is outlined in Algorithm 1, and its 



and V{a,x) ^ ^ 



convergence is established in [8]. 
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Fig. 2. <^{x,a*) vs. x 



Algorithm 1 Iterative Algorithm for Computing {a* , x* {a*)} 
States: 

Procedure: 

while \xk — Xk-i\ > e do 

cTfc-i = argmax{;7((7,a;fe_i) - Xk-iV{u,Xk-i)} 

0<cr<l 
U(ah-l,Xk-l) 



end while 



estimation accuracy plays an important role in the throughput 
performance: when a decreases from 1 to 0.1, the performance 
improves over 150% for p = 10. 

In Table I, we examine the convergence of the iterative 
algorithm I with a = 1. As expected, x{n) approaches to 
X* usually within a few iterations. 

TABLE I 

Convergence behavior of the iterative algorithm, a = 1. 



IV. Numerical Results 

In this section, we provide numerical examples to illustrate 
the above results. Unless otherwise specified, we assume that 
r, T, p, and M are chosen such that S = 0.1 and Ps = 
exp(-l). 
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Fig. 1. $(cr) vs. a, X = 0.1. 

Figure 1 depicts ^{x, a) as a function of the backoff ratio 
a. It can be seen that the average throughput is zero at both 
a = and cr = 1, and is maximized somewhere in between. 

Figure 2 depicts <I)(.t,(t*) as a function of x. Note that 
the optimal throughput x* is the solution to the fixed point 
equation x = ^{x,a*). It can be observed that x* is an 
increasing function of p for a given a, and is a decreasing 
function of a for a fixed p. It can also be seen that the 
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Table II compares the convergence behavior of the iterative 
algorithm with different error variance a and p — 1. When 
the error variance is large, the iterative algorithm needs more 
iterations to converge. Moreover, the backoff ratio a would 
decrease as a increases. This can be further observed in Fig. 
3. It indicates that when the estimation error is large, the 
transmitter would back off more to avoid channel outage. 

TABLE n 

Convergence behavior of the iterative algorithm, p = 1. 
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Table III illustrates the throughput gain, g = ^ ~^ , as a 
function of p, where x^ = $(0, a*) is the average throughput 
obtained by the schemes without using optimal scheduling. It 
can be seen that the throughput gain is more significant in the 
low SNR region, and is a decreasing function of p. 

In Table IV, we illustrate the throughput gain as a function 
of a. Note that a = 1/(1— /?)— 1 is an increasing function of (3. 
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Fig. 3. Backoff factor ct as a function of normalized error variance a 

TABLE III 
Throughput gain of DOS, a = 1. 



p 


0.5 


1 


2 


5 


10 


100 


X* 


0.254 


0.301 


0.336 


0.364 


0.374 


0.385 




0.185 


0.224 


0.254 


0.278 


0.288 


0.298 


9(p) 


37.3% 


34.3% 


32.3% 


30.9% 


29.8% 


29.2% 



As expected, when the normalized noise variance a increases, 
the optimal throughput x* decreases, as well as the x^. 
However, it is interesting to observe that the throughput gain 
increases instead. The rationale behind is that the performance 
of the schemes that do not use optimal scheduling "suffers" 
more than that of the distributed opportunistic scheduling in 
the presence of noisy channel estimation. 

TABLE rv 
Throughput gain of DOS, p = 0.5. 
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We also examine the performance of distributed oppor- 
tunistic scheduling as a function of the training time r. 
According to linear estimation theory, it has been shown that 
the error variance (3 and the training time r have the following 
relationship [8]: 

P = ^-r- (28) 

Using (28) in the simulation, we plot throughput performance 
of the distributed opportunistic scheduling as a function of 
the training time r in Fig. 4. It is clear that there exists an 
optimal training time which balances the tradeoff between 
better estimation accuracy and loss of transmission time. It 
can also be observed that when the average SNR p increases, 
the optimal training time decreases. 

V. Conclusion 
In this paper, we have generalized the study in [7] to 
investigate distributed opportunistic scheduhng under noisy 
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Fig. 4. Throughput vs. training time r 

channel estimation. For such cases, we have proposed that 
the transmission rate backs off on the estimated rate so as to 
reduce the outage probability. We have showed that the optimal 
scheduling policy has a threshold structure, but the threshold 
turns out to be a function of the variance of the estimation 
error, and is a functional of the backoff rate function. Since the 
optimal backoff is analytically intractable, we have proposed 
a suboptimal linear scheme that backs off on the estimated 
SNR (and hence the rate). The corresponding optimal backoff 
ratios and rate thresholds can be obtained using an iterative 
algorithm based on fractional maximization. Simulation results 
indicate that DOS still achieves significant throughput gain 
in the presence of noisy channel estimation, especially in 
the low SNR region. In addition, we have observed that 
the performance loss of DOS due to the imperfect channel 
estimation is less than that of the schemes that do not use 
channel-aware scheduling, indicating that the devised DOS is 
more robust against noisy channel estimation. 
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